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ABSTRACT 

Context. The extrasolar planet HAT-P-8b was thought to be one of the more inflated transiting hot Jupiters. 

Aims. By using new and existing photometric data, we computed precise estimates of the physical properties of the 

system. 

Methods. We present photometric observations comprising eleven light curves covering six transit events, obtained using 
five medium-class telescopes and telescope-defocussing technique. One transit was simultaneously obtained through four 
optical filters, and two transits were followed contemporaneously from two observatories. We modelled these and seven 
published datasets using the jktebop code. The physical parameters of the system were obtained from these results 
and from published spectroscopic measurements. In addition, we investigated the theoretically-predicted variation of 
the apparent planetary radius as a function of wavelength, covering the range 330-960 nm. 

Results. We find that HAT-P-8b has a significantly lower radius (1.321 ± 0.037 i?,iup) and mass (1.275 ± 0.053 A/jup) 
compared to previous estimates (1.50io o6 ^Jup 

and l.52to\l 

Afjup respectively). We also detect a radius variation in 
the optical bands that, when compared with synthetic spectra of the planet, may indicate the presence of a strong 
optical absorber, perhaps TiO and VO gases, near the terminator of HAT-P-8b. 

Conclusions. These new results imply that HAT-P-8b is not significantly inflated, and that its position in the planetary 
mass-radius diagram is congruent with those of many other transiting extrasolar planets. 
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1. Introduction 

' Transiting extrasolar planetary systems are of interest and 
I importance, because precise measurements of their phys- 
I ical properties can be achieved using spectroscopic and 
photometric observations. Atomic and molecular absorp- 
tion within the atmosphere of transiting extrasolar planets 
(TEPs) can also be investigated through transmission spec- 
troscopy (e.g. Swain et al. 2008; Sing et al. 2009; Fossati et 
al. 2010) and simultaneous multi-colour photometry (e.g. 
Ballester et al. 2007; Southworth ct al. 2012b) of the tran- 
sits. High-quality photometric observations not only en- 
able the measurement of the masses and radii of TEPs 
to accuracies of a few percent (e.g. Torres et al. 2008; 
Southworth 2009), but also the detection of transit anoma- 
lies due to stellar pulsations (Collier Cameron et al. 2010), 
tidal distortion (Li et al. 2010; Leconte et al. 2011), addi- 
tional bodies (moons, planets) (Kipping et al. 2009; Tusnski 
& Valio 2011), gravity darkening (Barnes 2009; Szabo et 
al. 2011) and star spots (Pont et al. 2007; Rabus et al. 
2009; Desert 2011). Besides the Rossiter-McLaughlin effect 
(Queloz et al. 2000; Gaudi et al. 2007), the photometric 
follow-up on consecutive/close nights of transits of plan- 



ets over parent-star starspots represents another fascinat- 
ing method (Sanchis-Ojeda & Winn 2011; Sanchis-Ojeda et 
al. 2011, 2012; Tregloan-Reed et al. 2012) to measure the 
sky-projected spin-orbit alignment. 

The increasing number of TEPs discovered every year is 
progressively revealing a remarkabel diversity. The improv- 
ing statistical weight of this sample is useful for establish- 
ing the correct theoretical framework of planet formation 
and evolution. Accurate estimates of the planet properties 
(mass, radius, orbital semi-major axis, etc.) are vital for 
this purpose, and photometric follow-up of known TEPs 
can dramatically improve our knowledge of the planet's 
characteristics. 

HAT-P-8 b is a transiting hot Jupiter found by the 
HATNet team (Latham et al. 2009), orbiting with a period 
of ~ 3.07 days around a star of spectral type F8 (Jones 
& Sleep 2000) or F5 (Bergfors et al. 2012). At the time 
of its discovery it was labelled as one of the most inflated 
transiting giant planets, with a measured mass and radius 
of Mb = 1.52i°;i^Mjup and i?b = 1.50t° i?jup, respec- 
tively. These values differ by 2~3f7 from the theoretical pre- 
dictions of Fortney et al. (2007). 
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The Rossiter-McLaughlin effect fias been detected in 
the HAT-P-8 system using radial velocity observations 
from the SOPHIE and FIES spectrographs. Simpson et 
al. (2011) found a sky-pro jcctcd orbital obliquity of A = 
-9.7ly:7 degrees ^nd Moutou et al. (2011) found A = 
— 171^^5 degrees; both values are consistent with alignment 
between the orbital axis of the planet and the rotational 
axis of the star. The two studies between them suggested 
lower values for Mb and R^, but neither calculated the 
physical properties of the system. 

Bergfors et al. (2012) have found a faint companion to 
the HAT-P-8 system using lucky im,aging observations with 
the AstraLux Norte instrument at the Calar Alto 2.2 m 
telescope. The companion, a likely M2-4 dwarf, is at an 
angular distance of 1.027 ± 0.011 arcsec and is fainter in 
the SDSS i' and z' passbands by Ai' = 7.34 ±0.10 mag and 
Az' = 6.68±0.07mag. The faintness of this star means that 
it has a negligible effect on optical observations of HAT-P-8. 

In this work we present eleven new follow-up light curves 
six transits in the HAT-P-8 system, obtained using five 1.2- 
2.5 m telescopes. We augment these data with previously 
published observations of seven transits, and measure the 
physical properties of the system. We find a substantially 
lower mass and radius for the planet, removing its outlier 
status and relegating it to a more well-populated part of 
the planetary mass-radius diagram. 

Our paper is structured as follows. In Sect. 2 we de- 
scribe the observations and data reduction. In Sect. 3 we 
analyse the data, and in Sect. 4 we obtain refined orbital 
ephemerides and physical properties of the HAT-P-8 sys- 
tem. In Sect. 5 we investigate the variation of the plan- 
etary radius as function of wavelength. Several anomalies 
detected in the light curves are discussed in Sect. 6, whereas 
in Sect. 7 we summarize the results and draw our conclu- 
sions. 



2. Observations and data reduction 

Six transits of HAT-P-8 b were monitored by five differ- 
ent telescopes between 2009 and 2012. Two transits were 
followed simultaneously by two of the telescopes, and one 
was simultaneously observed through four optical filters. 
Except the first, all our transit observations were performed 
with the telescope defocussing method, in order to min- 
imise the effect of Poisson, scintillation and flat-fielding 
noise (Southworth et al. 2009a) . We used autoguiding dur- 
ing all the observations, and the stars incurred a drift of 10 
pixels or less over each observing sequence. The correspond- 
ing night logs are reported in Tablet, and the differential 
photometry is tabulated in Table 2 and plotted in Figs. 1 
and 2. 

2.1. Kuiper 1.55 m telescope 

We observed one transit of HAT-P-8 in October 2009 using 
the University of Arizona's 1.55 m Kuiper telescope on Mt. 
Bigelow, Arizona. We used the Mont4k CCD, binned 3x3 
to 0.43'7pixel, for a total field of view (FOV) of 9.7' x 9.7', 
and an Arizona-/ filter^ . The observations were conducted 
with autoguidinp; and the telescope focussed. Due to the 

^ Tlic transmission curve for tiiis filter is shown at; 
james.as.arizona.edu / 'psmith/61inch /FILTERS /harris.jpg 



bright nature of the star we used 5 s exposure times, which 
resulted in an observing cadence of «30s. 

Systematic effects were minimised by autoguiding, keep- 
ing star wander to less than 5 pixels (2.15") over the course 
of the night. The resulting images were bias-subtracted, 
flat-fielded, and bad pixel-cleaned in the usual manner. 
Aperture photometry was performed using an IDL^ pipeline 
utilising the FIND and APER (Stetson 1987) tasks available 
in the NASA Astronomy User's Library''. The size of the 
aperture was chosen to minimise scatter in the data and was 
10 pixels (4.3") in radius. Several combinations of reference 
stars were considered, and the one which gave the lowest 
scatter in the final light curve was adopted. The 1083 orig- 
inal datapoints were binned to yield 216 final datapoints. 

2.2. Cassini 1.52m telescope 

One transit of HAT-P-8 was observed in November 2009, 
two in October 2011 and one in October 2012, using the 
1.52 m Cassini Telescope at the Astronomical Observatory 
of Bologna in Loiano (Italy). We have previously used this 
telescope several times to observe planetary transits (e.g. 
Southworth et al. 2012a), with the BFOSC (Bologna Faint 
Object Spectrograph & Camera) instrument operated in 
imaging mode. 

The CCD was used unbinned, giving a plate scale of 
0.58"/pixel, for a total FOV of 13' x 12.6', and the tele- 
scope was autoguided and defocussed. The first transit was 
observed through a Gunn r filter, the two 2011 transits 
through a Gunn i filter, and the last one through a Johnson 
/ filter. The first 2011 transit suffered from systematic noise 
due to a bad pixel in the aperture of the target star, so 
we did not use these data in our analysis. The transit of 
2012 was disturbed by clouds, which affected the photom- 
etry particularly at the end of the transit. We removed the 
points compromised by the clouds. 

The observations were analysed using the idl pipeline 
from Southworth et al. (2009a). The images were debiased 
and flat-fielded using standard methods, then subjected to 
aperture photometry using the APER task. Pointing varia- 
tions were followed by cross-correlating each image against 
a reference image. We chose the aperture sizes and com- 
parison stars which yielded the lowest scatter in the final 
differential-photometry light curve. The relative weights of 
the comparison stars were optimised simultaneously with 
fitting a second-order polynomial to the outside-transit ob- 
servations in order to normalise them to unit flux. 

2.3. Isaac Newton Telescope 

One transit of HAT-P-8 was monitored using the 2.5 m 
Isaac Newton Telescope (INT), La Palma (Spain), equipped 
with the Wide Field Camera (WFC) at prime focus. We 
used only one of the four CCDs, unbinned and with a plate 
scale of 0.33"/pixel, for a total FOV of 12.6' x 11.3'. The 
telescope was defocussed and autoguided, and the observa- 
tions were obtained through a Stromgren y flltcr. A few dat- 

^ The acronym IDL stands for Interactive Data 
Language and is a trademark of ITT Visual 
Information Solutions. For further details see 
http: //www. ittvis . com/ProductServices/IDL. aspx. 

The IDL Astronomy User's Library (ASTROLIB) is avail- 
able at http://idlastro.gsfc.nasa.gov/ 
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Fig. 1. Phased light curves of HAT-P-8 compared to the best jktebop fits using the quadratic LD law (left-hand panel). They 
are ordered according to central wavelength of the filter used. The residuals of the fits are plotted in the right-hand panel, offset 
to bring them into the same relative position as the corresponding best fit in the left-hand panel. 



apoints at the start of the observing sequence were rejected 
as they are affected by systematic noise (due probably to 
high airmass >2). The observations were reduced in the 
same way as those from the Cassini Telescope (Sect. 2.2). 

2.4. Calar Alto 1.23 m telescope 

Three transits of HAT-P-8 b, two in October 2011 and one 
in October 2012, were obtained using the 1.23 m telescope 
at the German- Spanish Calar Alto Observatory (CAHA) 
near Almerfa, Spain. Autoguiding was used. During the 
2011 observations, we used the 2kx2k SITE#2b otical 
CCD, which had^ a FOV of 16' x 16' and a pixel scale 
of 0.5" per pixel. We defocussed the telescope, read out 
only a small window in order to limit the dead time be- 
tween exposures, and observed through a Johnson R filter. 
The first transit was incomplete due to clouds, so was not 
included in our analysis. The 2012 transit was obtained 

* This CCD has been decommissioned. 



through a Johnson / filter with the new DLR-MKIII cam- 
era, which is equipped with an e2v CCD231-84-NIMO-BT 
DD sensor with 4kx4k pixels and a FOV of 21' x 21' at 0.3" 
per pixel. Unfortunately, the operator did not included ap- 
propriate reference stars in the field of view, and the scatter 
of the resulting light curve is higher than those of 2011. All 
the observations were reduced as for the Cassini Telescope 
(Sect. 2.2). 

2.5. Calar Alto 2.2 m telescope 

We observed one full transit of HAT-P-8 on the night of 
2012 August 26, using the 2.2 m telescope and BUSCA im- 
ager at CAHA. BUSCA is designed for simultaneous four- 
colour photometry: the light is split into four wavelength 
bands from UV to visual IR using three dichroics. In the 
four corresponding focal planes the same area of the sky is 
imaged onto 4kx4k 15/im pixel CCDs. For our observations 
we chose to have the SDSS u filter in the bluest arm and 
standard Calar Alto Gunn g, r and z filters in the other 
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Table 1. Details of the observations presented in this work. A'oba is the number of observations, Moon is the fractional illumination 
of the Moon at the midpoint of the transit, and toxp is the exposure time in seconds. The aperture sizes are the radii in pixels of 
the software apertures for the star, inner sky and outer sky, respectively. Scatter is the r.m.s. scatter of the data versus a fitted 
model, in mmag. Target mean PSF area in px^ is also reported for each dataset. Times and dates are in UT. /3 is the ratio between 
the noise levels due to Poisson noise and to combined Poisson and red noise. 



Telescope 


Date 


Start/End times 


Nobs 


^cxp 


Filter 


Airmass 


Moon 


Aperture 


Scatter 


PSF 




Kuipor 


2009 10 22 


04 


03 09 


10 


1083 


5 


Cousins I 


1.00 1.9 


25% 


10, 20, 25 


1.54 






Loiano 


2009 11 12 


16 


52 22 


09 


103 
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Gunn r 


1.02 1.01 1.55 


23% 
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1.15 


1963 


1.4 
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1.4 
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70 
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61% 
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0.83 


1662 


1.0 


CA 2.2m 


2012 08 26 


20 


31 04 


48 


160 


80 


sdss u 


1.52 1.00 -s- 1.55 


77% 


25, 55, 75 


2.37 


2827 


1.0 


CA 2.2m 


2012 08 26 


20 


31 04 


48 
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80 


Gunn g 


1.52 1.00 -i- 1.55 


77% 


30, 50, 70 


0.97 


2463 


1.2 
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Fig. 2. Phased BUSCA hght curves of HAT-P-8 compared to 
the best jktebop fits using the quadratic LD law (left-hand 
panel). They are ordered according to central wavelength of the 
filter used (sdss u, Gunn g, Gunn r, Gunn z). The residuals of 
the fits are plotted at the base of the figure, offset from zero. 



three arms. This choice led to a reduced field of view (from 
12' X 12' to a circle of 6' in diameter), but had the advan- 
tage of a much better throughput in grz compared to the 
default Stromgren filters. We defocussed BUSCA in such a 
way as to have as much signal as possible in the u band 



Table 2. Excerpts of the light curves of HAT-P-8: this table 
will be made available at the CDS. A portion is shown here for 
guidance regarding its form and content. 



Telescope 


Filter 


BJD (TDB) 


Diff. mag. 


Uncertainty 


Loiano 


r 


2455148.210781 


0.0009980 


0.0018456 


Loiano 


r 


2455148.214867 


0.0000050 


0.0017916 


Loiano 


i 


2455840.256418 


0.0000785 


0.0016988 


Loiano 


i 


2455840.258316 


0.0006256 


0.0016925 


CAHA 


R 


2455840.319832 


0.0007127 


0.0013972 


CAHA 


R 


2455840.329288 


-0.0019752 


0.0013798 


INT 


V 


2455437.412356 


0.0052976 


0.0009779 


INT 


y 


2455437.414575 


0.0056092 


0.0009998 


Kuiper 


I 


2455126.66906 


-0.99507 




Kuiper 


I 


2455126.66930 


-0.99959 





whilst remaining in the linear regime in the other pass- 
bands. The CCDs were binned 2 x 2 to shorten the read- 
out time. The autoguider was operated in-focus. The ob- 
servations were reduced in the same way as those from the 
Cassini Telescope (Sect. 2.2) and the resulting light curves 
are plotted in Fig. 2. 



3. Light curve analysis 

Sect. 2 introduced eleven transit light curves which were 
suitable for detailed analysis. These were each fitted using 
the JKTEBOP^ code (Southworth et al. 2004), which rep- 
resents the star and planet as biaxial spheroids. The pri- 
mary parameters of the fit were the orbital inclination, i, 
the sum and ratio of the fractional radii of the star and 
planet, ta + J'b and k = r\,/rA, defined as ta — Ra/o, and 
fh — Rh/o., and transit midpoint, Tq. The orbital semima- 
jor axis is a and, Ra and i?b are the absolute radii of the 
two celestial bodies. 

Once a fit was available for each dataset, we rescaled 
the errorbars to give a reduced of xl — 1- This step 
is necessary because the APER aperture photometry pro- 
cedure has a tendency to underestimate the measurement 
errors. Then, in order to take systematic noise into account, 
we inflated the errorbars further using the /3 approach (e.g. 
Gillon et al. 2006; Winn et al. 2008, 2009; Gibson et al. 
2008; Nikolov et al. 2012). We calculated /3 values for be- 



^ JKTEBOP is written in FORTRAN77 and the source code is 
available at http://www.astro.keele.ac.uk/~jkt/ 
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tween two and ten datapoints for each light cuurve, and 
adopted the largest (3 vahic. 

Limb darkening (LD) was accounted for using a 
quadratic law. The linear LD coefficient was fitted, whereas 
the non-linear one was fixed at a theoretically predicted 
value (Claret 2004), but perturbed by ±0.1 during the 
process of error estimating. The atmospheric parameters 
of the star assumed for deriving the limb- darkening co- 
efficients were: T^s = 6130, log.g = 4.15, [Fe/H]=+0.01, 
V^icro = 2 km s~^. Uncertainties in the fitted parameters 
from each solution were calculated in two ways: from 1000 
Monte Carlo simulations and with a residual-permutation 
algorithm (see Southworth 2008). The larger of the two pos- 
sible error bars was retained in each case. The light curves 
and their best-fitting models are shown in Figs. 1 and 2, 
whereas the parameters of each fit are reported in Table 3. 

We also attempted to fit all light curves simultaneously 
using the tap Gazak et al. (2012) code, but were unsuc- 
cessful. Analyses of individual light curves, however, yielded 
similar results and uncertainties as to those from JKTEBOP. 

3.1. Datasets taken from literature 

Latham et al. (2009) reported three z-band transits (1451 
points in total), only one of which was observed in its en- 
tirety, obtained with the 1.2 m telescope and KeplerCam 
at the F. L. Whipple Observatory, US. We converted 
the timestamps into orbital phase (using their ephemeris), 
sorted them and binned them into 154 points. 

Moutou et al. (2011) presented five datasets, two of 
them obtained by the Adagio Association with an 82 cm 
telescope {V band; 720 points), one from the 1.20m tele- 
scope at Observatoire de Haute-Provence {V band; 256 
points), one from a 35cm telescope at the Oversky obser- 
vatory (i? band; 569 points binned in 113 points), and one 
from a 32 cm telescope at the Quarryview Observatory {R 
band, 157 points). The two Adagio datasets were phased- 
binned into 144 points before analysis, to save computing 
time. The Oversky data were also binned, into 113 points. 

One transit was observed with the 36 cm Universidad 
dc Monterrey Observatory (UDEM) telescope (/ band; 639 
points). These data were reported by Todorov et al. (2012). 
Finally, one very good light curve was obtained by F. 
Hormuth using the CAHA 1.23 m telescope {R band; 176 
points). These data are available from TRESCA web site 
and were already studied by Simpson et al. (2011). 

All of the light curves in this Section were fitted using 
JKTEBOP in the same manner as our own data. 

3.2. Final photometric parameters 

All the light curves we fitted are shown in Figs. 1 and 2. 
The parameters of the fits are given in Table 3. In order 
to calculate the final photometric parameters we took the 
weighted mean of the individual values and uncertainties. 
This process highlighted the poor agreement between the 
different datasets. The xt of the parameters of the individ- 
ual light curves with respect to the final weighted means 
are 1.5 for ta, 2.4 for i and 2.7 for r^. The worst agreement 
is found for k, where xt = 6.9. Analagous situations have 
been found many times in the course of the Homogeneous 
SiAidies project (Southworth 2008, 2009, 2010, 2011, 2012), 
but the agreement is rarely this poor. 



Inspection of Fig. 1 gives a clue to this problem. The 
residuals of the fits to almost all of the available light 
curves, including our own, exhibit systematic deviations 
from zero. This correlated noise means the amount of in- 
formation available in the observational data is lower than 
suggested by the number of datapoints and their uncertain- 
ties. We have already increased our errorbars due to corre- 
lated noise using the /? approach (Sect. 3), but this was not 
enough to solve the problem. Fortunately, we have 18 sepa- 
rate light curves which have been analysed independently. 
Making the assumption that correlated noise in data taken 
on different nights using different telescopes is independent, 
we can account for it by inflating the errorbars on the final 
photometric; parameters by y^xE^ to give = 1-0 for each 
parameter. Because our results are based on 18 different 
datasets, we are confident that our analysis has yielded re- 
liable parameters. The light curve anomalies are discussed 
further in Sect. 5. 

3.3. Orbital period determination 

The transit time for each dataset was obtained using JK- 
TEBOP, and uncertainties were estimated by Monte Carlo 
simulations. In the determination of the orbital period of 
the HAT-P-S system, we also considered 23 timings mea- 
sured by amateur astronomers and available on the ETD® 
website (see Table 4). 

We excluded from the analysis the incomplete ETD light 
curves and those with a Data Quality index higher than 3. 
All timings were placed on BJD(TDB) time system. The re- 
sulting measurements of transit midpoints were fitted with 
a straight line to obtain a new orbital ephemeris: 

To = BJD(TDB)2 454 437.6742(14) + 3.0763458(24) E, (1) 

where E is the number of orbital cycles after the reference 

epoch (which we take to be the midpoint of the first transit 
observed by Latham et al. 2009) and qiiantities in brackets 
denote the uncertainty in the final digit of the preceding 
number. The fit has = 5.30, and the uncertainties given 
above have been increased to account for this. The large x^ 
indicates that the uncertainties in the various To measure- 
ments are too small, most probably due to the systematic 
differences between the light curves and their best fits as 
noted in Sect. 3.2. A plot of the residuals around the fit is 
shown in Fig. 3 and does not indicate any clear systematic 
deviation from the predicted transit times. We therefore 
take the conservative option of not interpreting the large 
X^ as a suggestion of transit timing variations. 

4. The physical properties of HAT-P-8 

Following Southworth (2009), we estimated the physical 
properties of the HAT-P-8 system from the photometric 
parameters measured in Sect. 3, published spectroscopic re- 
sults, and theoretical stellar models. The orbital eccentric- 
ity, velocity amplitude and metallicity of the star, (e = 0, 
Kp, ^ 153.1 ± 3.9, [Fe/H]= +0.01 ± 0.08) were taken from 
Latham et al. (2009), and for the parent-star effective tem- 
perature (Teff = 6130 ± 80) we used that measured by 
Kmitson et al. (2010). 

® The Exoplanet Transit Database (ETD) website can be 
found at http://var2.astro.cz/ETD 
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Table 3. Parameters of the jktebop fits to tfie ligfit curves of HAT-P-8. Tlie final parameters, given in bold, are the weighted 
means of the results for the 11 datasets. Results from theliterature are included at the base of the table for comparison. 



Source 


Filter 


+'"b 


k 


i° 


rA 




Loiano #1 


Gunn r 


0.177 ±0.013 


0.0921 ± 0.0034 


87.3 ± 2.3 


0.162 ±0.012 


0.0149 ±0.0014 


Loiano #2 


Gunn i 


0.1695 ±0.0030 


0.0943 ± 0.0022 


89.2 ± 1.0 


0.1460 ±0.0026 


0.01460 ± 0.00046 


Loiano #3 


Johnson / 


0.1674 ±0.0035 


0.0884 ± 0.0010 


88.5 ± 1.0 


0.1538 ±0.0032 


0.01360 ± 0.00034 


CA 1.23 m #1 


Johnson R 


0.218 ±0.012 


0.1006 ±0.0021 


82.7 ± 1.1 


0.198 ±0.010 


0.0199 ±0.0012 


CA 1.23 m #2 


Johnson _R 


0.1685 ±0.0061 


0.0888 ± 0.0032 


88.7 ± 1.6 


0.1548 ±0.0058 


0.01374 ± 0.00084 


CA 1.23 m #3 


Johnson / 


0.208 ±0.018 


0.1056 ±0.0042 


83.3 ± 1.6 


0.188 ±0.016 


0.0199 ±0.0018 


CA 2.2 m 


sdss u 


0.199 ±0.024 


0.0993 ± 0.0039 


84.1 ± 2.5 


0.181 ±0.021 


0.0180 ±0.0026 


CA 2.2 m 


Gunn g 


0.1720 ±0.0054 


0.0926 ± 0.0024 


89.9 ± 1.4 


0.157 ±0.050 


0.01458 ± 0.00057 


CA 2.2 m 


Gunn r 


0.1692 ±0.0051 


0.0959 ± 0.0018 


88.7 ± 1.4 


0.1544 ± 0.0046 


0.01480 ± 0.00059 


CA 2.2 m 


Gunn 2 


0.1695 ± 0.0059 


0.0917 ± 0.0016 


88.3 ± 1.5 


0.1552 ± 0.0053 


0.01424 ± 0.00058 


INT 


Stromgren y 


0.1699 ±0.0052 


0.0886 ± 0.0016 


87.9 ± 1.2 


0.1560 ±0.0047 


0.01383 ± 0.00056 


Kuipcr 


Cousins / 


0.175 ±0.012 


0.0871 ± 0.0026 


86.4 ± 1.6 


0.161 ±0.012 


0.0140 ±0.0010 


UDEM 


Cousins / 


0.220 ±0.032 


0.0836 ± 0.0039 


81.7 ±2.6 


0.203 ±0.030 


0.0170 ±0.0021 


Adagio 


Cousins V 


0.192 ±0.017 


0.0956 ± 0.0034 


85.3 ±2.1 


0.175 ±0.015 


0.0167 ±0.0018 


OHP 


Cousins V 


0.172 ±0.031 


0.0977 ± 0.0088 


87.4 ±3.9 


0.157 ±0.027 


0.0153 ±0.0031 


Quarryview 


Red filter 


0.169 ±0.011 


0.0901 ± 0.0024 


87.5 ± 2.1 


0.155 ±0.010 


0.0140 ±0.0011 


Oversky 


Sloan r' 


0.171 ±0.016 


0.0888 ± 0.0033 


86.9 ±2.6 


0.157 ±0.015 


0.0140 ±0.0016 


KeplerCam 


Sloan 2 


0.1781 ±0.0081 


0.0955 ± 0.0013 


86.9 ± 1.3 


0.1626 ±0.0073 


0.01553 ± 0.00084 


Final results 






0.09208 ±0.00049 


87.08 ±0.36 


0.1590 ±0.0014 


0.01468 ±0.00017 


Latham et al. (2009) 
Simpson et al. (2011) 




n 179C: + U.UUU94 

U.i (ZO_Q QQQ4g 

n 1700+0.00060 

U.i /00_Q QQQgg 


0.0953 ± 0.0009 
0.09135 ± 0.00089 


87.5±i;'^ 
87.80to;77 


U.10(0_Q QQ42 

iRS4+'''°°53 

U.iDJ4_g QQgg 


01501+"-"""'-*= 

0149'^+" """'*3 
U.Ui4yj_o 00066 



>1 
•a 



(a) 



<U -0.005 - 
Pi 
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OC) 
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Fig. 3. Panel (a): plot of the residuals of the timings of mid-transit of HAT-P-8 b versus a linear ephemeris. The timings in black 
are from this work, orange from Latham et al. (2009), green from Moutou et al. (2011), red from Todorov et al. (2012), and blue 
from Hormuth (TRESCA). The other timings obtained by amateur astronomers are plotted using open circles. Panels (b) and (c) 
are zooms in to the two best sampled regions. 



An initial estimate of the velocity amplitude of the 
planet (ifb) was iteratively refined by calculating the sys- 
tem properties using standard formulae (e.g. Hilditch 2001), 
and comparing the observed T^ft and ta with values of Toff 
and Ra/o. predicted by theoretical models for the calcu- 



lated mass of the star. This calculation was performed over 
a grid of ages and for five different sets of stellar models 
(see Southworth 2010). Statistical errors were propagated 
by a perturbation analysis. The resulting estimates of the 
physical properties are given in Table 5. The final set of 
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Table 4. Transit mid-times of HAT-P-8 and their residuals. TRESCA and AXA refer to the "TRansiting ExoplanetS and 
CAndidates" and "Amateur Exoplanet Archive" websites, respectively 



Time of minimum Cycle Residual 

BJD(TDB) -2400000 no. (JD) 



Reference 



Latham et al. (2009) 

Srdoc (AXA) 

Kucakova H. (TRESCA) 

Norby (AXA) 

Ayiomamitis (AXA) 

Srdoc (AXA) 

Kucakova H. (TRESCA) 

Trnka J. (TRESCA) 

Srdoc (AXA) 

Tieman B. (TRESCA) 

This work (Kuiper 155 cm) 

This work (Loiano 152 cm) 

Todorov et al. (2012) (Udem 36 cm) 

Moutou et al. (2011) (OHP 120 cm) 

This work (INT 250 cm) 

Moutou et al. (2011) (Oversky 35 cm) 

Hormuth F. (CA 123 cm) 

Ruocco N. (TRESCA) 

Marino G. (TRESCA) 

Moutou et al. (2011) (Hose 32cm) 

Hose K. (TRESCA) 

Dfeveny R., Kuchtak B. (TRESCA) 

Brat L. (TRESCA) 

Brat L. (TRESCA) 

Trnka J. (TRESCA) 

Zibar M. (TRESCA) 

Marek P. (TRESCA) 

Shadic S. (TRESCA) 

Brat L. (TRESCA) 

Trnka J. (TRESCA) 

This work (C A 123 cm) 

This work (Loiano 152 cm) 

Dvorak S. (TRESCA) 

This work (CA 220 cm BUSCA-w) 

This work (CA 220 cm BUSCA-g) 

This work (CA 220 cm BUSCA-r) 

This work (CA 220 cm BUSCA-z) 

Hose K. (TRESCA) 

This work (Loiano 152 cm) 

This work (CA 123 cm) 



54437.67461 ± 0.00044 0.00046 

55034.4824 ± 0.0031 194 -0.00284 

55034.4881 ± 0.0017 194 0.00286 

55046.7904 ± 0.0008 198 -0.00022 

55071.400 ± 0.002 206 -0.00138 

55071.4017 ± 0.0009 206 -0.00032 

55071.4056 ± 0.0012 206 0.00422 

55071.4132 ± 0.0013 206 0.01182 

55074.4752 ± 0.0013 207 -0.00253 

55123.7021 ± 0.0018 223 0.00284 

55126.77446 ± 0.00079 224 -0.00115 

55148.3140 ± 0.0019 231 0.00397 

55409.7992 ± 0.0012 316 -0.00022 

55434.4108 ± 0.0023 324 0.00061 

55437.48610 ± 0.00043 325 -0.00043 

55437.4834 ± 0.0011 325 -0.00313 

55437.48434 ± 0.00088 325 -0.00219 

55437.4895 ± 0.0010 325 0.00297 

55440.5646 ± 0.0012 326 0.00171 

55449.79074 ± 0.00099 329 -0.00118 

55452.86716 ± 0.00079 330 -0.00110 

55797.4152 ± 0.0018 442 -0.00379 

55797.4186 ± 0.0015 442 -0.00039 

55800.4916 ± 0.0013 443 -0.00374 

55800.4942 ± 0.0011 443 -0.00114 

55800.4968 ± 0.0025 443 0.00146 

55800.4996 ± 0.0010 443 0.00426 

55812.80021 ± 0.00085 447 -0.00051 

55834.3372 ± 0.0019 454 0.00206 

55837.4136 ± 0.0025 455 0.00211 

55840.48811 ± 0.00063 456 0.00028 

55840.48845 ± 0.00049 456 0.00062 

55886.6345 ± 0.0021 471 0.00148 

56166.5776 ± 0.0011 562 -0.00289 

56166.58178 ± 0.00057 562 0.00129 

56166.58162 ± 0.00051 562 0.00113 

56166.57924 ± 0.00054 562 -0.00124 

56175.81123 ± 0.00059 565 0.00170 

56203.49662 ± 0.00031 574 -0.00002 

56203.49682 ± 0.00083 574 0.00018 



physical properties was obtained by taking the unweighted 
mean of the five sets of values found from the different stel- 
lar models, which also allowed us to obtain an estimate of 
the systematic errors inherent in the use of stellar theory. 
The results of this process are given in Table 6. 

Finally, following the method delineated by Enoch et al. 
(2010) and improved by Southworth (2011), we used em- 
pirical measurements of stars in detached eclipsing binary 
(dEB) systems to calibrate the parent star of the HAT-P-8 
system. This allowed us to measure the physical proper- 
ties of the system without using stellar models, thus avoid- 
ing the dependence on stellar theory. These results are also 
given in Table 5. 

Table 6 also contains a comparison between our own re- 
sults and those found by Latham et al. (2009). We find 
smaller masses and radii for both the planet and the host 
star. This is surprising because the values of the fractional 
radius of the star, which furnishes the vital constraint on 
the stellar density, arc very similar between the two studies 
(we find VA = 0.1590 ± 0.0014 versus ta = 0.1575t^;^j^89 for 
Latham et al. 2009). One small difference is that we adopted 
Teff = 6130 ± 80K (Knutson et al. 2010) as opposed to the 
value of 6200 ± 80K used by Latham et al. (2009). Our re- 
vised physical properties move HAT-P-8 b into a region of 



parameter space which is comparatively well-supplied with 
transiting planets. Fig. 4 shows the change in position in 
the planet mass-radius plot. 

We contacted D. Latham in order to check the veracity 
of the final values of the parameters that they reported in 
their discovery paper. D. Latham confirmed that some of 
these results are indeed correct, and has kindly supplied re- 
vised values. These are from a re-analysis carried out by J. 
Hartman on a slight updated datasct, and arc given in the 
final column of Table 6. They correspond to a smaller plan- 
etary mass and radius (Mb = 1.38 Mjup, = 1.40i?jup) 
than given in Latham et al. (2009), in agreement with our 
own findings. 



5. Variation of planetary radius with wavelength 

One of the factors that plays a principal role in deter- 
mining the atmospheric properties of hot- Jupiter planets 
is the amount of stellar flux incident on the planet's sur- 
face. Variations in this irradiation cause different planets 
to have different atmospheric chemical mixing ratios and 
atmospheric opacities. This could potentially lead to divide 
into classes of hot Jupitcrs. An initial suggestion was that 
of Fortney et al. (2008) to distinguished pM- and pL-class 
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Table 5. Derived physical properties of the HAT-P-8 planetary system. 





This work 


This work 


This work 


This work 


This work 


This work 




(dbB constraint) 


(Claret models) 


(Y^ models) 


(Teramo models) 


(VKss models) 


(UsLP models) 


Kb (kms-^) 


159.1 ±3.9 


154.8 ±2.3 


156.7 ± 1.9 


153.4 ± 1.7 


154.2 ± 2.6 


155.4 ±2.13 


Ma (Mr^) 


1.292 ± 0.094 


1.189 ± 0.054 


1.234 ± 0.044 


1.158 ± 0.037 


1.175 ± 0.061 


1.202 ± 0.053 


Ra (Rq) 


1.521 ±0.041 


1.480 ±0.032 


1.470 ± 0.026 


1.467 ± 0.023 


1.474 ±0.031 


1.485 ±0.026 


logflA (cgs) 


4.1852 ±0.0137 


4.1733 ±0.0088 


4.1954 ±0.0103 


4.1694 ±0.0105 


4.1715 ±0.0106 


4.1748 ±0.0115 


Mb (Mjup) 


1.345 ±0.071 


1.273 ±0.048 


1.305 ± 0.041 


1.251 ±0.037 


1.263 ±0.053 


1.282 ±0.047 


Rb (Rjup) 


1.357 ±0.043 


1.320 ±0.034 


1.337 ± 0.032 


1.308 ±0.030 


1.315 ±0.037 


1.325 ±0.034 


Ph (Pjup) 


0.503 ±0.034 


0.517 ± 0.034 


0.511 ±0.033 


0.522 ±0.034 


0.519 ±0.035 


0.516 ±0.034 


e 


0.0691 ± 0.0026 


0.0711 ±0.0023 


0.0702 ± 0.0022 


0.0717 ±0.0022 


0.0714 ±0.0025 


0.0708 ± 0.0023 


a (AU) 


0.04510 ±0.00109 


0.04387 ± 0.00066 


0.04442 ± 0.00053 


0.04348 ± 0.00047 


0.04370 ± 0.00075 


0.04403 ± 0.00065 


Age 




4 7+1.4 

^- ' -0.5 


o c:+0.4 


4 7+1-0 
^■'-0.5 


4 2+l'5 


4 2+°'^ 



Table 6. Final physical properties of the HAT-P-8 system, compared with results from the literature. Where two errorbars are 
given, the first refers to the statistical uncertainties and the second to the systematic errors. 





This work (final) 


Latham et al. (2009) 


Moutou et al. (2011) 


Latham (private comm.) 


Ma (Mq) 


1.192 ±0.061 ±0.043 


1.28 ±0.04 




1 r,o+U.U4 
^•^°-0.06 


Ra (Rq) 


1.475 ± 0.032 ± 0.010 


1 co+0.08 
J^-OO_o.06 




1.57 ±0.07 


logffA (cgs) 


4.177 ±0.011 ±0.019 


4.15 ±0.03 




4.15 ±0.03 


PA (p©) 


0.371 ±0.013 ±0.018 






0.46 ± 0.05 


Mb (Mjup) 


1.275 ± 0.053 ± 0.030 


1 rn + 0.18 


1.34 ±0.05 


1.38 ±0.05 


Rh (Rjup) 


1.321 ±0.037 ±0.016 


1 ■S0+°°^ 

l.UU_g Qg 




1.40 ± 0.08 


gb (ms-^) 


18.11 ±0.82 


16.98 ±1.17 




1.38 ±0.05 


Pb (pjup) 


0.517 ±0.035 ±0.006 


0.568 ± 0048 




0.62 ± 0.09 


(K) 


1713 ± 24 ± 13 


1700 ± 35 




1771 ± 39 


e 


0.0710 ± 0.0025 ± 0.0008 


0.061 ± 0.003 




0.069 ± 0.004 


a (AU) 


0.04390 ± 0.00075 ± 0.00052 


0.0487 ± 0.0026 


0.0449 ± 0.0007 


U.Uy:OU_o.oo07 


Age (Gyr) 


A q+1.5+0.8 

^•■^-1.4 -0.5 


3.4 ± 1.0 




q q+0.7 
■^"J-O.S 



in their atmospheres. By observing a planetary transit at 
different wavelengths, it is possible to detect a variation 
in the value of the radius measured as a function of the 
wavelength, and probe chemistry and wavelength depen- 
dent opacity at the planet's terminator. This dependence 
on the wavelength can be therefore used to probe the at- 
mospheric composition of TEPs, being a complementary 
method to the observations of secondary eclipses. 

As an additional possibility offered by the BUSCA 
data, we made an attempt to investigate possible varia- 
tions of the radius of HAT-P-8 b in different optical pass- 
bands. Receiving from its parent star an incident flux of 
2.22±0.20x 10^ ergs^^ cm^^ (Latham, private communica- 
tion), HAT-P-8 b should belong to the pM class of planets. 
The theoretical models of Fortney et al. (2010) predict that 
its radius should be few percent lower at 350-400 nm and 
800-950 nm versus 500-750 nm. Following the strategy used 
by Southworth et al. (2012b), we fitted the BUSCA light 
curves with all parameters fixed to the final values reported 
in Table 4, with the exception of k and the LD coefficients. 
The errors were estimated by a residual-permutation algo- 
rithm. The results are exhibited in Fig. 5, where the points 
show the data, the vertical bars represent the relative er- 
rors in the measurements and the horizontal bars show the 
full widths at half maximum transmission of the passbands 
used. As is apparent from Fig. 2, the BUSCA-z light curve 
shows some structure in its residuals from phase to 0.02, 
which suggests that the z-band data suffer correlated noise. 
The result coming from this band should be therefore con- 
sidered with caution. The optical region not covered by the 
BUSCA data was investigated by using the i-band light 
curve from Loiano, which is the best one that we have at 
this band (see Fig. 1). 



2.5 




OOLi , , , , 1 

0.0 0.5 1.0 1.5 2.0 2.5 

Planet mass (Mjup) 

Fig. 4. Plot of the masses and radii of the known TEPs. The 
orange symbols denote values from the Homogeneous Studies 
project and the blue symbols results for the other known TEPs. 
HAT-P-8 b is shown in red (Latham et al. 2009) and green (this 
work). Dotted lines show where density is (from bottom to top) 
1.0, 0.5, 0.25 and 0.1 pjup. 



planets, depending on the presence of strong absorbers such 
as gaseous titanium oxide (TiO) and vanadium oxide (VO) 
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A variation of rb along the five passbands is clearly visi- 
ble in Fig. 5. We investigated this variation with the help of 
model atmospheres. We used detailed non-gray atmosphere 
codes to model the temperature structure and transmis- 
sion spectrum of the planet. We computed ID model atmo- 
sphere of HAT-P-8b, using the atmosphere code described 
in Fortney et al. (2005, 2008). Pressure-temperature pro- 
files are derived that either include or exclude the opacity 
of TiO and VO molecules. The fully non-gray model uses 
the chemical equilibrium abundances of Lodders & Fegley 
(2002) and the opacity database described in Freedman et 
al. (2008). The atmospheric pressure-temperature profiles 
simulate planet-wide average conditions or day-side aver- 
age conditions. We computed the transmission spectrum of 
the models using the methods described in Fortney et al. 
(2010) and Shabram et al. (2011). 

Excluding TiO and VO from the opacity calculation, in 
the upper panel of Fig. 5 we compare the transmissions spec- 
trum of the ID planet-wide average profile (in blue) with 
a warmer day-side average model (in green) , to experimen- 
tal data. In comparison, the bottom panel of Fig. 5 shows 
analogous calculations in red and yellow respectively, which 
include both the TiO and VO species. Coloured boxes indi- 
cate the predicted values for these models integrated over 
the bandpasses of the observations. Models without TiO 
and VO have optical transmission spectra that are dom- 
inated by Rayleigh scattering in the blue, and pressure- 
broadened neutral atomic lines of Na an K at 589 nm and 
770 nm, respectively. Models with TiO and VO absorption 
show significant optical absorption (a much large transit 
radius), that broadly peaks around 700 nm, with a sharp 
falloff in the blue, and a shallower falloff in the red. 

Comparing the panels of Fig. 5, it is readily apparent 
that the model that gives the best match to the data is 
the one with TiO and VO opacity in the atmosphere of 
HAT-P-8 b. The increased optical radius is somewhat larger 
than the model prediction. Taken at face value, the obser- 
vations are in general agreement with the Fortney et al. 
(2008) 's hot- Jupiter classification based on stellar irradia- 
tion. We suggest that HAT-P-8 b should be an important 
target for foUowup studies to confirm or refute these sugges- 
tive observations. The clear detection of TiO/VO, or other 
strong optical absorbers, would be an important step in 
characterizing hot Jupiter atmospheres, as such absorbers 
are thought to cause temperature inversions in these at- 
mospheres (e.g. Fortney et al. 2008; Burrows et al. 2008). 
Previously, Desert et al. (2008) suggested a detection of 
TiO and VO at strongly depleted levels in the red optical 
spectrum of HD 209845 b, using Hubble STIS. 

6. Light-curve anomalies 

While one of our follow-up light curves, obtained at the 
INT, displays a regular transit shape, the others all show 
anomalies of a similar shape. They display an asymmetry 
with respect the line of minimum transit time. It is im- 
portant to clarify if the anomalies that we detected have 
an astrophysical nature or are of instrumental or environ- 
mental origin. One way to check this is to have independent 
measures of the same transit event, obtained using multiple 
telescopes located at different sites. In this way, if both the 
telescopes reported the same anomaly, it less likely that 
they are caused by instrumental or Earth- atmosphere ef- 
fects. This already happened inadvertently in the case of 



0.0154 
0.0152 
0.0150 
0.0148 
0.0146 
0.0144 
0.0142 
0.0140 

0.0154 
0.0152 
0.0150 
0.0148 
0.0146 
0.0144 
0.0142 
0.0140 
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Fig. 5. Variation of the fractional planetary radius rb — Rb/a 
with wavelength. The points shown in the plot are from the 
Calar Alto 2.2 m (u, g, r and z bands) and Loiano (i-band) 
telescopes. The vertical bars represent the errors in the mea- 
surements and the horizontal bars show the full widths at half 
maximum transmission of the passbands used. The experimen- 
tal points are compared with four models. These use planet- 
wide average pressure-temperature profiles (in blue and red) or 
warmer day-side average profiles (in green and yellow). Synthetic 
spectra in the top panel do not include TiO and VO opacity, 
while spectra in the bottom do, based on equilibrium chemistry. 
Coloured squares represent band-averaged model radii over the 
bandpasses of the observations. 



several follow-up observations of WASP-33 carried out by 
amateur astronomers and reported in Kovacs et al. (2012). 
It also occurred by chance for the transit observations of 
HAT-P-8 on the night of 2010/08/28, which was observed 
by ourselves at the INT and contemporaneously at CAHA 
by F. Hormuth. Fig. 6 shows the two light curves in the 
same plot; we note that both have a regular transit shape 
but disagree over the transit depth. Some fraction of this 
disagreement is due to the different LD characteristics, as 
the INT data were obtained with a bluer filter than the 
CAHA data. 

In 2011, we aimed to observe transits in the HAT-P- 
8 system from two different observatories and this goal 
was successfully achieved on the night of 2011/10/05 when 
a transit was simultaneously observed using the 1.52 m 
Cassini and CAHA 1.23 m telescopes. The datasets show 
partial but not complete agreement about the transit shape 
anomalies, as well as slight differences due to the different 
LD in the i and R passbands. Details of the two light curves 
are displayed in the lower panel of Fig. 7, while the upper 
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1.005 



1.000 



-a 
'3 



0.995 

a 

0.990 



0.985 



55 437.3 55437.4 55437.5 55437.6 

BJD(TDB) - 2400000 

Fig. 6. An example of the consistency of follow-up observations 
made of the same transit on the night of 2010/08/28 by different 
telescopes. This is an example of a "regular" transit. Lighter yel- 
low dots are for the data taken at the INT (this work), whereas 
the darker brown ones for those taken by F. Hormuth (ETD) at 
the CAHA 1.23 m. The agreement between the two datasets is 
not very good. The difference in depth between the two datasets 
is partly due to different filters used, Stromgren y and Johnson 
R, respectively. The vertical line represents the expected transit 
minimum time. 



panel shows another transit (2009, Loiano) which presents 
a similar anomaly. The anomalies cannot be removed by 
choosing different comparison stars for the differential pho- 
tometry, but are not completely consistent between differ- 
ent datasets for the same transit. We have therefore treated 
them as correlated noise in our analysis (see Sect. 3). 

If the recurrent anomalies that we noticed in the HAT- 
P-8 light curves are not of instrumental origin, is there a 
reasonable astrophysical explanation for these features? We 
can easily exclude a gravity darkening effect (e.g. Barnes 
2009) because the rotation rate of the star is low (Latham et 
al. 2009; Moutou et al. 2011). The presence of a moon is an 
unlikely possibility, and in any case is difficult to constrain 
with so few observations. 

Another possibility is that the anomalies are due to sin- 
gle or a belt of stellar spots on the photosphere of the par- 
ent star. Knutson et al. (2010) estimated that HAT-P-8 
has an activity index of log(i?^^) = —4.985, which indi- 
cates that the star has a moderate chromospheric activity. 
This is related to the star's magnetic structure and there- 
fore to the presence of photospheric features, such as spots, 
that modulate luminosity. It is thus possible that, dur- 
ing the transit, dark starspots are occulted by the planet. 
Similar features were indeed already observed in several 
cases, e.g. HD 189733 (Pont et al. 2007), TrES-1 (Rabus 
et al. 2009), Keplcr-17 (Desert 2011), HAT-P-11 (Sanchis- 
Ojeda & Winn 2011), and WASP-19 (Tregloan-Reed et al. 
2012). One argument against this is the relatively high Tcs 
of the star, 6130 ± 80 K, which makes starspots less likely 
to occur. Moreover, spots would only affect small parts of 
the light curve, while here both light curves show trends for 
the full duration of the transit. 

Another possible explanation is that the faint M-dwarf 
companion of HAT-P-S is a flare star, which could emit 
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Fig. 7. Upper panel: an example of a discrepant transit. The 
green points are for the data taken at the Cassini telescope (this 
work) through a Gunn r filter. Lower panel: an example of the 
consistency of followup observations made on the same date of 
2011/10/05 by different observers. Lighter blue dots are for the 
data taken at the Cassini telescope (this work) , whereas the red 
ones for those taken at CAHA 1.23 m (this work). Some of the 
difference in depth between the two datasets is due to different 
filters used, Gunn i and Johnson R, respectively. The vertical 
line represents the expected transit minimum time. 



giant flares bright enough to significantly affect the light 
curves. However, the magnitude difference between the two 
objects (Ai' = 7.34 ± 0.10, Bergfors et al. 2012), is such 
that this star would have to generate a superfiare in order 
to make a noticeable difference. One more explanation for 
this behaviour is that a significant fraction of the compar- 
ison stars display intrinsic variability which is insufficient 
to detect in individual stars but, when combined, is enough 
to modify the transit shape in the IIAT-P-8 system. 

Having exhausted all possibilities, we conclude that the 
apparent transit shape distortion shown in Fig. 7 could be 
caused by differential color extinction, pathological variabil- 
ity in the comparison stars, or other low-frequency noise of 
atmospheric or astrophysical origin. Further investigations 
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of these hypotheses would, inevitably, require substantial 
further data and is beyond the scope of the current work. 

7. Summary and Conclusions 

We have reported observations of six transits of the HAT- 
P-8 system performed using five different medium-class 
telescopes, for a total of eleven new light curves. All but 
one of these transits were obtained using the defocussed- 
photometry technique, achieving a photometric precision 
of 0.47-1.7mmag per observation. Four of them were ob- 
served on the same nights from two different telescopes. 
In one of these two nights, both the light curves show an 
anomaly which is probably caused by systematic noise of 
atmospheric or astrophysical origin. Another transit was si- 
multaneously observed through four optical passbands in a 
wavelength range between 330 and 960 nm. 

We modelled our new and seven published datasets us- 
ing the JKTEBOP code, and used the results to determine 
the physical properties of the planet and its host star. 
Compared to the discovery paper (Latham et al. 2009), we 
find a significantly smaller radius and mass for HAT-P-8b. 
We obtain i?b = 1.321 ± 0.037 i?jup versus 1.50t^;^^ i?jup, 
and Mb = 1.275 ± 0.053 Mjup versus 1.521°;^^ Mjup. The 
theoretical radius calculated by Fortney et al. (2007) for 
a core-free planet at age 4.5 Gyr and distance 0.045 AU 
is 1.107-1.108 iijup for a planet of mass in the range 1.0- 
1.46 Mjup, which is still significantly smaller than the radius 
we find. 

Latham et al. (2009) found that HAT-P-8 b was a com- 
paratively highly inflated planet. Instead, our results place 
it firmly in a well-populated part of the mass-radius dia- 
gram, removing its outlier status. Whilst the existing tran- 
sit light curves of the HAT-P-8 system all show some sys- 
tematic deviations from the best fits found using simple 
geometric models, the large number of available datasets 
means that our overall results are reliable. HAT-P-8 is an- 
other system where extensive follow-up photometry has 
been necessary to determine robust physical properties for 
a planetary system. 

Finally, thanks to the ability of BUSCA to measure stel- 
lar flux simultaneously through different fllters, covering a 
quite large range of optical window, we probed the compo- 
sition of the atmosphere of HAT-P-8 b by investigated how 
vary its radius in these wavelengths. In fact, the presence 
of strong optical absorbers in the atmosphere of the planet 
should produce larger transit radius at optical wavelengths 
than in the near UV or near infrared. We measured a vari- 
ation of the radius of HAT-P-8 b along five passbands, cor- 
responding to a wavelength coverage of 330 — 960 nm. This 
result was then theoretically investigated by using several 
synthetic spectra based on isothermal model atmospheres 
in chemical equilibrium. The comparison between the mod- 
els and the experimental points suggests the presence of 
molecular gas that strongly absorbs in the optical, poten- 
tially composed of TiO and VO gases, in the atmosphere 
of HAT-P-8 b. 
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